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The spatial and temporal evolution of parametrically excited geodesic acoustic mode 
(GAM) initial pulse is investigated both analytically and numerically. Our results show 
that the nonlinearly excited GAM propagates at a group velocity which is, typically, much 
larger than that due to finite ion Larmor radius as predicted by the linear theory. The 
nonlinear dispersion relation of GAM driven by a finite amplitude drift wave pump is also 
derived, showing a nonlinear frequency increment of GAM. Further implications of these 
findings for interpreting experimental observations are also discussed. 


I. INTRODUCTION 


Geodesic Acoustic Modes (GAM) [1, |2j] are finite-frequency components of zonal structures 


(ZS) [3, IJ] unique in toroidal plasmas, which could be spontaneously excited by microscopic drift 
wave (DW) type turbulence 5|, including drift Alfven waves (DAW); and in turn, are capable 
of scattering DW/DAW into stable short radial-wavelength domain Q, 0-0]. Therefore, GAM 
may regulate the turbulence intensity and the associated wave-induced transports llO, ll|. The 
excitation of GAM by DWs can be described by parametric decay instability fill fl3l]. where the 
DWpurnp (cuo,ko) decays into a DW lower sideband (w,s,k,s) and a GAM (wG,k<y), and selection 
rules of frequency and wavenumber matching conditions are satisfied. 

It is known that GAM has a finite linear group velocity due to finite ion Larmor radius (FILR) 
effects, and this linear group velocity is typically radially outward, consistent with GAM continuum 
due to radial temperature profiles. For GAM with k r pi <C 1, its linear group velocity is Vq = 


u 2 = lo 2 g { 1 + C G k 2 p 2 ) [2]. Here, 


log = 


C G io G k r p 2 , with the linear dispersion relation given as 
\/T e /Ti + 7/Avti/Ro is the local GAM frequency in the fluid limit, vu = \/2 Ti/nrii the thermal ion 
velocity, Rq is the major radius of the torus, k r is the radial wavevector, pi = rriiCVti/(eBo ) is the 
ion Larmor radius, and C G is an order unity coefficient. This linear group velocity of GAM has 
been discussed in several works [2], [ 3 ], and is shown to have important consequences on the 

nonlinear excitation of GAM by DWs and change the absolute/convective nature of the parametric 


instability [8], 16j]. Radial propagation of GAM has been observed in several experiments [17H19I]. 


and qualitative agreement between the experimental results and linear theory has been obtained 
[19l, 120], However, in-depth analysis of the experimentally obtained dispersion relation leads to the 
conclusion that, even though a quadratic dependence of GAM frequency on its radial wavevector, 
qualitatively consistent with linear theory of KGAM [jjj, is indeed obtained, the coefficient for 
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FILR effects is much larger than that predicted by linear theory [14, [ 21 ]. This discrepancy has 
also been found in numerical simulations 15j], where the measured radial propagation velocity of 


the DW driven GAM is used to determine the coefficients of FILR effect, and is found to be much 
larger than unity. Up to now, there is no first-principle-theory-based explanation of this “enhanced 
FILR effect”, although a general theoretical framework exists to formulate it with all necessary 
physics ingredients [8|. In fact, we will show in this work, that this discrepancy could be due to 
nonlinear effects. 

We note that, GAM is an n = 0/m ~ 0 mode, with m = ±1 sidebands and possibly higher order 
ones (depending on the perpendicular wavelength Q]), such that it is generally not driven unstable 
by expansion free energy of the plasma. Here, m and n are, respectively, the poloidal and toroidal 
mode numbers in the Fourier mode structure representation adopting straight field line toroidal 


flux coordinates [22J]. Thus, GAM, in general, could be observed when it is nonlinearly driven 
by ambient turbulence, and in this case, the spatial-temporal evolution of GAM is significantly 


affected by the DW nonlinear drive. It has been pointed out in Ref. [8j that the nonlinearly driven 
GAM propagates at a much larger nonlinearly-coupled group velocity in the presence of DW. As a 
result, the propagation of GAM and experimental observations should also be interpreted taking 
nonlinear effects into account. In this work, we shall further study the spatial-temporal evolution 
of the nonlinearly coupled DW-GAM system more in details in order to analyze its implications to 
experimental observations. The rest of the paper is organized as follows: the theoretical model used 
in this work is presented in Sec. QJ1 which is then solved in Sec. IHIl for the spationtemporal evolution 
of the coupled GAM-DW system. The possible applications of our theory to interpretation of 
experimental observations and numerical results are also presented. Finally, a brief summary is 
given in Sec. IIVI 


II. THEORETICAL MODEL 


The equations describing the nonlinear interactions between GAM and DW are derived using 
gyrokinetic theory [3, Q] . Assuming that the DW is constituted by a constant-amplitude pump 
wave and a lower sideband with a much smaller amplitude due to GAM modulation, the normalized 
coupled nonlinear equations describing GAM excitation by DW are then given as equations (9) 
and (10) of Ref. 

(d t + 7 s + iup - iu* - iC d u*p*dty A s = T^S, (1) 

[dt(dt + 2y g) + uq — CccoQPidr^ S = iLOG^odlAs- (2) 

Here, S = d r Ac/oi is related with the GAM electric held with a = ^(aiWpTe/Tj) 1 / 2 and at = 


1 + 5P±/(eno5cj)p ) an order unity coefficient [23|, <5Pi is the perturbed perpendicular pressure 


due to the DW pump scalar potential 8<j)p in the k±pi <C 1 limit j3[. Meanwhile, Ac, Ap 
and A$ are, respectively, the radial envelopes of GAM, DW pump and lower sideband, Tq = 




3 


{aiTi/upTeY^ckg^Ap/B is the normalized pump wave amplitude. Furthermore, 75 and 7 c are 
the Landau damping rates of DW sideband and GAM, wp is the pump DW frequency and w* is 
the diamagnetic drift frequency. The kinetic term in equation (jT|) , i.e., the term proportional to 
Cd, comes from finite radial envelope variation due to the coupling between nei ghb oring poloidal 
harmonics. The expression for Cd can be derived from equation (19) of Ref. 2J, and one has 


Cd ~ 0(e / {n 2 q' 2 pj)) with q being the safety factor, q' = dq/dr its radial derivative and e = t/Rq. 
On the other hand, the kinetic term in equation ([ 2 ]); i.e., the term proportional to Cq, comes 
from FILR of GAM. Thus, Cq ~ 0(1) and its detailed expression can be obtained from equation 
(9) of Ref. [21[. Other notations are standard. We note that, the governing equations (P) and 
(|2|) are derived from quasi-neutrality condition assuming both GAM and DW are predominantly 
electrostatic perturbations. Electrons respond adiabaticly to kn / 0 perturbations, i.e., DW and 
m/0 poloidal sidebands of GAM; while ion responses are solved assuming q ^$> 1, k±pi -C 1 and 
|w 0 1 ~ |w*| for DW. We note that, even though turbulence usually refers to a broad spectrum of 
nonlinearly interacting DWs, in the present analysis we have considered the nonlinear interactions 
of GAM with a single-n DW in order to elucidate the nonlinear effects on the radial propagation 
of GAM due to interaction with finite-amplitude DWs. Since for each DW with toroidal mode 
number n, the interactions with the corresponding GAM is coherent, we may expect that in the 
presence of DW turbulence consisting of multiple-n modes, the net nonlinear effects would be an 
appropriate sum/integral of the nonlinear effects of individual -n mode sconsidered here. System 
nonuniformities in equations (HJ and (J 2 ]) , which may affect qualitatively the convective/absolute 
nature of the parametric process as shown in Ref. @, are also ignored here in order to focus on the 
radial propagation of the parametrically excited GAM pulse. 


III. SPATIALTEMPORAL EVOLUTION OF THE COUPLED DW-GAM SYSTEM 


Equations (jT]) and ([ 2 ]) can be solved using two-spatial two-temporal scales expansion of S and 
As, i.e., As = As(t, () exp(ikor — iuot) and S = S(r, () exp(ikor — iuj^t) such that dt = — iwo + d T 
and d r = iko + d £, with r and ( denoting the slow temporal and spatial variations. In order 
to delineate the physics of GAM propagation, we can assume that the system is well above the 
excitation threshold [g]. Thus, we can ignore 75 and 7 g in equations (jTJ) and ([ 2 ]); and the coupled 
nonlinear equations reduce to 


(9 r + U s 5 c )i s = r * 0 S, 

(d T + Vcd c )S = ^r 0 (k 2 0 -2ik 0 d c )A s . 


( 3 ) 

( 4 ) 


Here, Vs = 2CdOJ*pfko and Vq = Cc^GPiko are, respectively, the linear group velocities of DW 
sideband and GAM. We note that Vs and Vq have the same sign for typical tokamak parameters 


|, such that the excitation of GAM by DWs is a convective amplification process, ignoring 
system nonuniformities ful l^ . Furthermore, in deriving equations © and ©, the following 
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frequency and wavenumber matching conditions for resonant decay are applied 

-w 0 + wp-w* + C d u]*klpj = 0, 

— LUq + UJq + Cc^G^OPi = O’ 


from where (wo, &o) can be solved for. 

Moving into the wave frame by taking £ = Q—V c t , with V c = (Vs + Vg)/ 2, the coupled nonlinear 
equations, ([3]) and (|4|). can be combined to yield the following equation describing the nonlinear 
spatialtemporal evolution of the parametrically excited GAM 

(o 2 - v 0 2 dl) i = l -k 2 v 2 i - ikovld^i. ( 5 ) 

Here, Vq = (Vs — Vq)/ 2. Letting S = exp(i/3£)A(£, r), with j3 = /coLg/(2V r 0 2 ), equation ([5]) reduces 
to 


( d 2 r - V 2 d 2 ) A = Qfcgrg + /3fc 0 r 2 - /3 2 H 0 2 ) A = f A. (6) 

Physically, /3 can be interpreted as nonlinear modification to the GAM wave vector, which also 
affects the GAM frequency as shown below. Equation ([6]) can be solved, and yields the following 
unstable solution 


A = 


J dk / exp f j exp ik^ + \]fj 2 - kjVtfr 


^/irAk 0 J- 

This solution corresponds to the initial condition 


Aq 


( 7 ) 


A = Aq exp — 


Akie 


at r = 0. We note that this is the typical wave packet initial structure for parametrically excited 
GAM, with a spectrum width Ako- Assuming |14<9,e| -C |<9 r |, i.e., convective damping due to FILR 


effects are higher order corrections to the temporal growth |2|, l27j], the general solution, equation 


CD. can then be reduced to the following time asymptotic solution: 

Aq 


A = 


AkoX-, 


■ exp TjT 


e_ 

4A2 


( 8 ) 


Here, A 2 = (1/A/cq + Vqt/( 2fj)), and it corresponds to GAM initial pulse broadening in time. The 
time asymptotic solution of GAM electric field is then 


S = 


Ai 


o 


Ako\ T 


exp ?r + i/3(( ~ V c t) - (( - V c t ) 2 ^j . 


(9) 


One then readily has from equation ([9]) that the nonlinearly excited GAM is characterized by 
a nonlinear radial wavevector 


k NL = k 0 - id c lnS = k 0 (l + T 2 /(2H 0 2 )) , 


(10) 
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i.e., the wavevector increases with pump DW amplitude, and is larger than that predicted from 
frequency/wavenumber matching conditions. 

The real frequency of the excited GAM can also be obtained from equation Q 


ojnl — wo T ^r ^ ^ — wq T 


koT 2 0 V c 

W 0 2 


( 11 ) 


o;o(fco) can be solved from the matching conditions, which can then be substituted into equation 
m , and yield: 


, k 0 T 2 0 V c 1^, ,22 

unl = w G + 9 ^ 2 + -C G UGk 0 Pi 


'0 
^ 2 i 


— W G H--h 


CG^GPik 2 NL 


2V 2 


( 12 ) 


2{l + T 2 /(2V 2 )r 

This is the nonlinear dispersion relation of the parametrically excited GAM. We note that, both 
Vo and V c are proportional to ko, and thus, the nonlinear frequency shift due to the modulation of 
DW, koT^Vd (2 Vq), is independent of k o. Thus, finite amplitude DW will increase the frequency of 
the nonlinearly driven GAM. The frequency increment, can be expressed as (ed(j)/T) 2 (L n / p^) 2 from 
our theory, which indicates an order of unity frequency increment for typical tokanrak parameters. 
This may explain the existence of the higher frequency branch of the “dual-GAM” observed in 
HT-7 tokanrak |19|], which oscillates at a frequency much higher than other branch with the usual 


GAM frequency (The frequencies of the two co-existing “dual - GAMs” are respectively 12 and 21 


kHz in shot 113901 |19(). Another finding of the HT-7 experiment is the coefficient of FILR effect 
is 0( 10 2 ) larger than that predicted by linear theory 28j. On the other hand, equation (1121) shows 
that the coefficient for kinetic dispersiveness is, in fact, decreased by a factor (1 + Tq/ (2Vq )) 2 . 
The reason why experimental analysis found an “increased” coefficient is that, in the analysis of 
experimental data, one employed the linear dispersion relation of GAM and used the expression 
(uJobs ~ Giioc) / (ui oc k 2 bs pj) to determine the coefficient C G 28[. Here, the subscript “o6s” denotes 
experimental observation, and “Zoc” denotes local continuum frequency of GAM. As we have shown 
in equation (1 12 [) . u co 0 b s — ui 0 ” contains the kinetic dispersiveness as well as the order one nonlinear 
frequency increment koT^Vj {2 Vq)\ which, thus, can lead to an over-estimation of the coefficient 
Cq L 28]. The effective coefficient obtained in this way is, Cq ~ (eS4>/T) 2 (L n /pi) 2 / (k^p 2 ) ~ 
O(10 2 ) for typical tokamak parameters; which is significantly larger than that predicted by linear 
theory. Our nonlinear theory, thus, provides a possible explanation of experimental observations. It 
can also be used to explain the O(10 2 ) increase of the FILR coefficient from numerical simulations 


m 


The coupled GAM and DW sideband wavepacket, propagates at a nonlinear group velocity 
V c = (Vs + Vg)/ 2, which is much larger than the linear group velocity of GAM due to |Rs| | V G \ 
(|o;p| ~ |w*| |w G | for resonant decay). Thus, to interpret the propagation of GAM nonlinearly 


excited by DW turbulences including DAW, linear theory of KGAM p, 21] is not adequate, and 
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FIG. 1: Nonlinear wavenumber k r v.s. pump amplitude 

r 0 


one must instead, apply nonlinear theory. We note also that, while both the real frequency and 
wavevector of the excited GAM depend on the amplitude of the pump DW, the nonlinear group 
velocity is determined by ko from matching conditions, and is independent of the pump amplitude. 
Thus, for the comparison of experimentally observations with analytical theory, the nonlinear group 
velocity may be a better candidate. 

The coupled nonlinear GAM and DW sideband equations, equations (JT]) and (O, are solved 
numerically. Here, we fix Cd = Cq = 1, wg = 0 . 1 , wp = w* = 1 , and study the coupled nonlinear 
equations by varying To- The dependence of the nonlinear wavenumber k r on pump amplitude 
is given in Fig. [TJ where the dots are the wavenumbers from numerical solution, the diamonds 
are the wavenumbers obtained from equation (HOI) ; and the solid curve is obtained from matching 
condition. For the parameters we have here, the wavevector solved from matching conditions is 
k 0 = 0.32. We may see from Fig. [T] that our nonlinear theory fits well with the numerical results; 
and it reduces to ko as To approaches 0. The comparison of the numerically measured nonlinear 
group velocity with our theory, is presented in Fig. 0 where the dots are numerical results and 
the diamonds are obtained from V c = (Vs + Vq)/ 2, and Vs and Vq are defined with ko- We 
note that, for the parameters we used in numerical solution, Vg = 0.64, Vq = 0.032 and V c = 
(Vs + Vq)/ 2 = 0.34 2> Vg- Very good agreement between numerical results and analytical theory 
(< 3% discrepency) are obtained here, suggesting that experimentally observed radial propagation 
of GAM must be understood using nonlinear theory. 

The nonlinear frequency of GAM is given in Fig. [3j where the dots are numerical results 
and the diamonds represents Unl from equation m■ Note that, for the parameters we use 
here, ojq = 0.105, and the nonlinear frequency from numerical solution increases with pump DW 
amplitude as predicted by our theory. 
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FIG. 2: Nonlinear group velocity v g v.s. pump 
amplitude Tq 



FIG. 3: Nonlinear GAM frequency a i NL v.s. pump 
amplitude Tq 


IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, equations describing the spatialtemporal evolution of parametrically excited GAM 
initial pulse are studied both analytically and numerically. It is found that the parametrically 
excited GAM propagates at a nonlinear group velocity, which is the mean of the linear group 
velocities of GAM and DW, and is much larger than that predicted by linear theory of kinetic 
GAM. The wavevector of the excited GAM has a quadratic dependence on the amplitude of the 
constant-amplitude pump DW. On the other hand, the nonlinear group velocity is independent of 
the pump DW amplitude; suggesting it as a good candidate for the comparison between experiments 
and analytical theory. Our nonlinear theory, further shows that there is a nonlinear upshift in 
the GAM frequency. Implications of the present theoretical findings to the HT-7 experimental 
observations are also important. Our results demonstrate that one must include nonlinear effects 
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in order to properly analyze numerical simulations and/or experimental observations of GAM. 

We note that, while the ambient turbulence in experiments consists of a whole spectrum of 
nonlinearly interacting DWs, we have, in the present analysis, considered the modification of the 
GAM dispersion relation due to a single-n DW with finite amplitude. This can be justified, since 
DW interactions with ZS, e.g. GAMs, have two components: a coherent part due to the interaction 
with the self-generated ZS, and a random contribution due to interaction with ZS produced by other 
incoherent components of the fluctuation spectrum [2a]. In both cases, the coupling coefficient is 
roportional to DW intensity and, therefore, we focus here on the coherent GAM-DW interaction 
3 ]. Effects of system nonuniformities, which are shown to play important roles on the nonlinear 
interactions between GAM and DWs are also ignored here. This will limit the applications of our 
nonlinear theory. To properly interprete global numerical simulations and/or experimental results, 
more in-depth investigations taking into account system nonunifomities will be needed. 
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